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FOREWORD 


This report describes a portion of the results obtained on NASA Grant 
NSG 3044. This work was done under subcontract to the University of 
Illinois, Urbana, with Prof. S.S. Wang as the Principal Investigator. I'ho 
prime grantee was the Massachusetts Institute of Technology, with Prof. 

F.J. McGarry as the Principal Investigator and Dr. J.F. Mandell as a 
major participant. The NASA ~ LeRC Project Manager was Dr. C.G. Chamis, 

Efforts in this project are primarily directed towards the development 
of finite element analyses for the study of flaw growth and fracture of 
fiber composites, This report presents a method of analysis for adhesive 
or Interlaminar cracks which propagate in the interface, rather than co- 
hesively in the adhesive or interlayer. The latter case was treated in an 
interim report NASA CR-135248. 


ABSTRACT 


A study on the elastic behavior o£ Interface cracks In adhesively 
bonded lap shear joints Is presented. The problem is Investigated by 
using a recently developed method of analysis based on conservation laws 
In elasticity for nonhomogeneous solids and fundamental relationships In 
fracture mechanics of dissimilar materials. The formulation leads to a 
pair of linear algebraic equations In mixed-mode stress Intensity factors. 
Singular crack-tip stress Intensity solutions are determined directly by 
Information extracted from the far field. Stress Intensity factors and 
associated energy release rates are obtained for various cases of Interest, 
Fundamental nature of the Interfaclal flaw behavior In lap-shear adhesive 
joints Is examined in detail. 
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1. INTRODUCTION 


While adhesively bonded joints have been employed in many secondary 
structural components, the effective use of adhesive joining technology 
in primary load-bearing structural members Is still In Its Infancy* Ap- 
plications of adhesive joints will not reach their full potential until 
certain critical problems are solved. The major factors responsible for 
this situation are the complex failure modes and mechanics presented in 
the joints. Because of many geometric, material and fabrication variables 
involved, the fracture problem of adhesive joints is extremely complicated. 
Current knowledge of the joint failure behavior is yet very limited. Ana- 
lyses and design criteria for flaw resistance of adhesively bonded joints 
are seriously deficient. 

Interfacial cracking, also called debonding, la one of the most fre- 
quently encountered modes of failure. It frequently occurs at geometric 
boundaries such as edges and re-entry corners due to inherent stress con- 
centrations, or results from faulty joining in fabrication such as in- 
complete wetting between adherend and adhesive. Debonding is also found 
in adhesive joints subjected to combined high temperature and moisture 
absorption. Progressive reduction of stiffness, exposure of the interior 
of the composite to environmental attack, and final disintegration of the 
structure are frequently observed to occur by the presence and growth of 
interfacial cracks. The lap-shear joint shown in Fig. 1 is considered In 
this study, because It Is one of the most widely used Joint configuration 
In structural applications. Understanding the fundamental nature of Inter- 
faclal cracks is of utmost Importance In the reliability and safety design 
of the adhesive joints for primary load-bearing structures. In this paper 
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analytical study on the interfacial crack behavior In the adhesively bonded 
lap-shear Joints Is presented. 

The Interface crack problem was apparently first studied by Williams 
[1]. Elasticity theories for cracks of this nature were attempted by many 
Investigators, and solutions characteristically involve oscillatory stress 
singularity [2-7]. A major concern for the oscillation of near-fleld 
stresses Is that they lead to physical absurdity of crack surface Inter- 
penetration or overlapping, as pointed out by England [2] and Malyshev, et 
al. [8]. This implies that the solutions for the Interface crack problem 
are physically Inadmissible. The unsatisfactory aspects of the oscillatory 
stress singularity were discussed by several researchers, for example, 

T/lefs. [2,8-10]. An alternate model of a closed crack tip with frlctlon- 
less contact was first proposed by Comnlnou [9] to correct these effects. 
The stress singularity of the model is different from that of the oscilla- 
tory one with a finite normal tensile traction ahead of the crack tip. The 
crack tip In a nominal tensile field has an extremely small contact area 
In comparison to the size of the crack. Recently, Achenbach, et al. [10] 
and Keer, et al. [11] proposed another models which Include the crack face 
closure, relative slip conditions at the Interface, and different defini- 
tions of stress Intensity factors. In view of the complexities aforemen- 
tioned, the Interfaclal crack problem in adhesively bonded joints is ob- 
viously rather difficult. The complex structural geometry, the presence 
of the multiphase material system, and the unknown multlaxlal stress state 
acting on the crack Introduce additional complications. The very local 
nature of the Interface crack singularity as shown in [7,12] Is noted to 
Introduce the thickness of the thin adhesive layer as a characteristic 
dimension In the problem. Any simplifications, which fall to consider the 
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crucial role of the thin adhesive layer, would lead to a severe drawback 
In studying the full nature of the adhesive joint fracture # The importance 
and complexities of interface debonding in adhesive joints have long been 
recognised [13-15] . However, research progress on this kind of real-life 
defect has been relatively slow. 

In this paper, a study on the inter facial crack behavior in adhesively 
bonded lap-shear joints is presented. By using the recently developed 
conservation laws of elasticity for nonhomogeneous solids [16-20] snd 
fundamental relationships in fracture mechanics of Interface cracks, the 
problem is reduced to a pair of linear algebraic equations, and stress 
intensity solutions can be determined directly by Information extracted 
from the far field. This feature makes the current approach particularly 
suitable and attractive. Solutions are obtained for adhesively bonded 
lap-shear joints with various material systems and geometric configurations 
of the crack and the joint. The method of analysis is of practical use in 
the design and analysis of adhesive joint fracture!. The fundamental nature 
and unique features of the interface flaw behavior in adhesively bonded 
joints are revealed. 


2, FORMULATION 


The Eshelby-Rice conservation law for a horoogeneous solid in a plane 
elasticity problem [16-19] has recently been extended to a solid composed 
of two dissimilar materials [20] with the following form: 


J.(s} 

1 


^ "1 - “jk "k ^ 


C[W]6i2 - 


where W is the strain energy density; the stress tensor; , the dis- 

placement vector; n^, the outward unit normal of an arbitrary closed contour 
S which encloses a portion of the continuum, and S, is a portion of the inter- 
face bounded by S as shown in Fig. 2. The bracket [ ] in Eq. 1 denotes the 
jump of a function across the interface C 2 * 0, i.e,, 


[W] - W(X 3 ^, 0 +) - W(x^,0-), 


( 2 ) 


(3) 


Continuity conditions of displacements and interlaminar stresses across 
X 2 * 0 require 

[Uj ] * 0 and [o^ 2 ^ “ 


The component of the conservation integral in Eq. 1 can be simplified to 


J^iS} 


' "1 - “jk "k 

Js 


(5) 


which is the same as the analogous result for a single phase material. 

For a solid containing an interface crack between two dissimilar materials 
as shown in Fig. 3, in Eq. 5 along a path T has the standard J-integral form 
in homogeneous fracture mechanics, i.e.. 
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( 6 ) 


(WdX2 

r 



ds), 


:e r ia an arbitrary path aurrounding the crack tip, provided that the 
± surfacea are free from traction and tt.e interface la a atraight line. 

:h Independent integrals (aimilar to Rice’s J^integral of course with appro- 
kte modifications) for elastic media with spatially varying moduli were noted 
Ltkinson [21] through an energy-momentum tensor forraaliam.] The J is shown 
relate to the energy release rate G of an interface crack in a usual manner, l.e., 

J » G - Gj + Gjj., (7) 


re Gj. and G^^ are the energy release rates associated with the mode I and 
! II stress intensity factors. The interface crack-tip stress intensity 
:ors, Kj. and are defined in a complex form as 



2^ a®" Urn (z) , 

Z-XJ 


» 


( 8 ) 


:e is a complex potential in the well-known Kolosov-Muskhelishvill 

nulation [3,22], and $ is a bimaterial constant given in Ref. [3]. 

It is noted chat and Kj.j, defined in Eq. 8 for an interfacial crack 

lissimilar media are different from those for a crack in a homogeneous 
Ld. Thus, the and may not possess the usual significance and 
jical interpretation as in the cohesive (or homogeneous) fracture. While 
overall stress Intensity factor [3,23] may be used to express the maximum 
Litude of stresses and to correlate crack extension, accurate description 
:he oscillatory singular stress field in the neighborhood of the crack tip 
jlres detailed knowledge of individual stress intensity factors, 
luse of the complexity of the problem, numerical procedures are necessary 
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obtaining accuratt solutions for Intarfacs crack probXoms. For examplSt 
and I'lar [24] used a hybrld’-s tress finite element method to achieve this, 
and Stern [23] employed a contour integral method based on Betti's 
;>rocal work theorem and, recently, Sroelaer [25] used crack">flank dis- 
ament data provided by numerical solutions to yield the interface crack- 
stress intensities. Each of these numerical approaches requires its own 
itatlonal scheme to handle the problem and haa given satisfactory results, 
lis section, an alternate and convenient method of analysis is proposed for 
mnining the inter facial and in bonded joints. 

Using the stress field established by Rice et al. [3], one can readily 
that the J«lntegral in Eq. 7 is related to the stress intensity factors 
an interfacial crack by 


J 



1 



(Kf + Kfj) 


(9) 


a plane strain case. The coefficients in Eq. 9 are replaced by 1/[4 m^ 

v^) ] for a plane stress condition. The J-integral alon& does not provide 
aate information for determining individual values of and for an 
rently, mixed-mode interfacial crack. Further development of Eq. 9 based 
le introduction of known auxiliary solutions for the crack problem can 
/e this difficulty. 

Consider two independent equilibrium states with associated field variables 
ted by superscripts 1 and 2 for the elastically deformed blmaterial 
. Superposition of the two equilibrium states leads to another equilibrium 
a, 0. The J-integral for the superimposed state can be shown to have the 


j(0) .,(1) ^ j(2) ^„a.2)_ 


( 10 ) 
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c« «noeh«r coni^rvAtlon lnt«gr«X with th« 0orm 


« 


( 1 , 2 ) . 


f 


^,(1,2) _ 


au<2) 31,(1) 

t ( 1 ) ^V- , M) S 
h +^i 3x 


dt 


( 11 ) 


w 


( 1 . 2 ) 


in Eq. 11 is the mutuml potenclal energy density of che bineteriel 


ft defined by 


" °ljkl l,j \,i “l,J “k,* 


( 12 ) 


illing the J-K relacionship in Eq. 9 for the superioposed stace and com- 
ing it with Eq. 10 » one can obtain 




(13) 


re 


a - y T~ (1 - V . ) * 
iil ‘♦^i ^ 


(14) 


M-integrala in Eqs . 11 and 13 deals with inceraccion terms only, and are 
d directly in solving the interface crack problem. The M-integral is 
arly related to the details of the stresses and deformation at the crack 
(i.e., and in Eq. 13), but yet may be evaluated in the far field 
e., the integral in Eq. 11), where such a calculation can be carried out 
ti greater accuracy and convenience than near the crack tip. It is also 
ed that, while ICj. and characterize the controversial near-field oscil- 
ory singular stresses, the energy release rate G, and, perhaps, and G^^, 
quantities well defined and can be evaluated conveniently mathematically 
physically. 

^Thc M-integral used here and elsewhere in this paper is defined by Eq. 11, 
should not be confused with those in Refs. [17} and [13] . 
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3. SOLUTION PROCEDURE 


Equation XX together with Eq. 13 provides sufficient information for 
:ermining the stress intensity solutions for a mixed-mode interface crack 
>blem» when known auxiliary solutions arc introduced* Denote the first 
ciliary solution by a superscript 2a for a crack between two dissimilar 
:erials subjected to mode-I loading only with 

- 1 and - 0 * (15) 


lation 13 can be simplified to 

jj(l,2a) ^ 2a (16) 

sre has the same form as Eq. 11 with the subscript 2 being replaced 

2a > T^^^ and u^^^ in Eq. 11 can be determined along a properly selected 
t^'.|iie&‘?ion path V in the far field by any convenient method such as the 
lunonly used finite element analysis. For a plane crack problem with the 
iding of Eq. 15, T^^^^ and u^^^^ are derivable by the Kolosov-Muskhelishvili 
emulation, and the auxiliary solution has the following form; 




f^^^(x,n r, 


9; V^, P2» ^^2^^ 


(17) 


l:h 


T,(2a) 

i 


n 

ij j ’ 


(18) 


- -p= t, 8, 

let forms of the functions f^^^ and g^^^ may be found in [3,4]. 

The second auxiliary solution, denoted by the superscript 2b, for the 
itaterial solid under a pure-mode II loading such that 
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K^.O «.d kW-1. (20) 

can be found by the same formulation. The conaervatlon Integral ia, there- 
fore , expreaaed by 

. 2a (21) 

where haa a form similar to that of Eq. 11. 

(2a) (2a) 

It Is Important to note that the auxiliary solutions, u^ ' and 

independent of the particular boundary-value problem 
being posed. Therefore, they may be determined Independently by any con- 
venient analytical method once for all. In solving for and the 

Integrals, and mxist be evaluated accurately and explicitly. 

For a given crack geometry, loading condition and bimaterial constant, this 
can be achieved by Integrating Eqs. 16 and 21 along a properly selected 
contour in the far field so as to avoid crack-tip complications. In 
conjunction with the auxiliary solutions determined, a numerical method, 
using a conventional finite-element approach, is currently employed to 
calculate T^^^ and u^^^ . The M-integrals are then formally evaluated by 
using the second-order Gaussian quadrature along a contour F passing through 
Gaussian stations of each element (Fig. 4). 

'1 
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4, NUMERICAL RESULTS AND DISCUSSION 


The solution procedure aforementioned is programmed for studying general 
two-dimensional Interface crack problems. Evaluation of the conservation 
integrals is conducted in conjunction with a conventional finite element 
method by using eight-node isoparametric elements. Solutions are obtained 
for lap-shear joints with various adherends, adhesives, and geometric con- 
figurations. Accuracy and convergence of the results are affected by several 
unusual features of the problem and the method of analysis due to the singular 
nature of the crack and inherent approximation involved in the numerical 
evaluation of the conservation integrals. Assessments of solution accuracy 
are made by examining relevant problems for which unquestionably correct 
and exact solutions are available in the literature. Excellent agreement 
is obtained between current results and reference solutions. Details of the 
study are reported elsewhere [26]. Current results have an accuracy within 
approximately three percent deviation from reference solutions based on the 
optimum mesh and the integration contour presently selected. The primary 
objectives of this section are to determine stress intensity solutions, 
which characterize local stresses and deformation, and to examine effects 
of material and geometric variables on the crack- tip response. Of particular 
interest are the complex failure modes and energy release rates associated 
with the interface crack. 

4.1 Interfacial crack and Joint Configuration 

The lap-shear joint considered in Fig. 1 is composed of two high- 
stiffness and high-strength adherends bonded by a thin adhesive layer. The 
upper and lower adherends and the adhesive layer are assumed to have uniform 
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thicknesses of t^ and respectively. The overlap region has a 

dimension 1 . A crack of length CL is located along the interface of the 
upper adherend and the adhesive near the traction-free edge. Except for 
the interfacial crack» perfect bonding is assumed everywhere. The two 
adherends are made of materials with elastic constants and 

and the thin adhesive layer has properties of E^ and ^2 • The adhe- 
sive and adherends are assumed to be linear, elastic. Studies of related 
cases such as the center-of-bond cohesive crack and the eccentric crack 
problems were reported elsewhere [27-28]. 

4.2 Effects of Adhesive Properties on Interfacial Crack-Tip Response 
Geometric variables in the problem are given as the following; 

6 ■ 45" , L ■ 0.5 in. , 

tj^ “ t^ ■ 0.05 in., t 2 ■ 0.005 in., (22) 

a * 2,5 t 2 ■ 0.0125 in. 

The upper and lower adherends are made of the same aluminum with elastic 
constants, E^^ ■ E^ ■ 10 x 10® psi, and " 0.33 . Effects of 

different kinds of adhesive on the crack-tip response are examined by 
considering various values of E 2 • For the purpose of generality, adhes- 
ives with a broad range of modulus values covering three decades on a 
logarithmic scale are studied. The cases with typical 

for aluminum/ epoxy systems. Higher E^^/E 2 values correspond to the joints 
with less rigid adhesives or subjected to a "hot and wet" environment. 

Crack-tip stress Intensity solutions and associated energy release rates 
are obtained as functions of the modulus ratio shown in. Figs. 5 and 6. 

Failure modes in the joints are observed clearly. The interfacial crack 
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oxperlances o mixed-mode fractute even under uniaxial loading. Both Kj 
and Kjj decrease rapidly as the adhesive modulus decreases, and the 

Increases with The mode~I stress 

intensity factor is higher than that of the raode-II in the entire range 
studied, except for the cases of ~ suggests that the doml~ 

nant mode of failure is opening-mode. In the case of E^/E 2 •“ 1 , the crack 
is cohesive, and the nature of the crack-tip singularity changes. The values 
of Kj and Kj.j in this case are found to be consistent with the results 
in [27], The total energy release rate and the energy release rate of in- 
dividual fracture mode are determined also. For given adhe rends , the 
reduction of E 2 increases the coefficients of and in Eq. 9, 

which influencq the G value significantly. The is found to increase 


difference between and 


rapidly with ^II* other hand, decreases as a less stiff 

adhesive is used. The rapid increase of G^. and G corresponds to a decrease 
of Kjj. and in the Interface crack problem — a unique phenomenon not 

observed in the homogeneous crack problem in general. In the case of very 
large adhesive may become incompressible with Poisson's ratio 

approaching 0.5. The incompressibility of the soft adhesive is not considered 
in the present study. 


4.3 Effects of Dissimilar Adherends 

Effects of dissimilar adherends on the interface crack behavior are 
studied for lap-shear joints with the same geometry and crack length used 
in the previous section. The adhesive layer and the lower adherend 
have the following elastic properties: E 2 ■ 0.5 x 10® psi, V 2 ■ 0*35 and 

Ej ■ 10 X 10® psi, ■ 0.33. The problem is solved by considering the 
joints with various ranging from E^/E 2 ■ 1 to 1000. The results 
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are given in Figs. 7 and 8, in which , Kj.j, G^, Gj.j. and G are related 
to E-/E- in aemi-logarithnilc plots. Adherends with dissimilar properties 

X M 

are found to have significant effects on the crack-tip response. The 
opening-mode stress intensity factor increases appreciably but the shearing- 
mode stress intensity remains relatively constant as the stiffness of the 
upper adherend approaches that of the lower one, l.e,, ^ 

E^/E 2 exceeds 20, an opposite trend of changes of Kj. and are 

observed; Kj. remains relatively unchanged and Increases gradually. 

The ratio of in the joint with dissimilar adherends is always 

smaller than, or at most equal to, that of a joint with the same upper and 
lower adherends. Effects of dissimilar adherends on energy release rates 
are shown in Fig. 8. The G^ is approximately one order of magnitude 
larger than G^^ in the entire range of ^ total energy 

release rate increases rapidly initially, and then remains relatively constant 
with the change of E^/E 2 . As the stiffness of the upper adherend exceeds 
that of the lower one, crack resistance of the joint becomes less sensitive 
to the dissimilarity of the two adherends. But a higher fracture resistance 
is obtainable as the upper adherend becomes less stiff. 

4.4 Effects of Adhesive-Layer Thickness 

Effects of the adhesive- layer thickness on fracture of the joint have 
long been recognized. Bascom, et al. , for example, [13,14] showed experi- 
mentally that fracture energy release rates are related closely to the 
bond thickness. Extensive analytical studies on debonded composite lami- 
nates were conducted by Erdogan [5,6,29]. Several unique features of the 
geometric variable on the composite failure behavior were revealed. Wang, 
et al. [27] studied center-of-bond cohesive cracks and reported fundamental 
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aracterlstlcs o£ thickness effects on the adhesive joint fracture. In 
U portion of the study, lap-shear Joints with a geometric configuration 
d crack length Identical to the one used in previous sections are considered, 
a joints are made of identical aluminum adherends and epoxy adhesive with 
Etstic constants: ■ 10 x 10^ psi, " ^3 " 0.33, and E 2 “ 0.5 x 

^ psi, " 0.35, respectively. The adhesive thicknesses, ranging from 
/40 to tj^/5 (i.e., from 0,00125 in. to 0,01 in,), are considered. 

Relatlonshlpfj among Kj., G^, G and the normalized bond thick- 

3 s, ^ 2 ^^! **1 “ *"3 " shown in Figs. 9 and 10. The Kj. 

found to be insensitive to the change of the adhesive-layer thickness, 
t increases rapidly with t 2 * Quantitatively, K^. is always larger than 
j. in all cases studied. This is particularly true for the joint with a 
inner adhesive, in which has a value about one order of magnitude higher 
an that of in general. The K^. clearly dominates the fracture of the 
int with a thinner bond, but the shearing mode failure becomes increasing- 
important as the adhesive thickness increases. Associated energy release 
tes are given in Fig. 10 for joints with various t 2 /tj^'s. The decrease 
the adhesive-layer thickness leads to an increase of fracture resistance 
the joints (i.e., a reduction in the total energy release rate and the 
ximum cleavage stress) . The change of fracture resistance with the bond- 
ne thickness in a lap-shear joint is consistent with the results reported 
Ref. [29]. However, variations of interfacial crack-tip stress intensity 
ctors and associated energy release rates are different from the solutions 
tained for center-of-bond cohesive crack problems, in which K’s and G's 
e almost independent of the adhesive-layer thickness. This situation is 
pected as stress singularities are different in the two cases. 
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4.5 Effects of Crack Length 

Effects of the Interface crack langth are of significant Interest 
because of the particular Joint configuration and the loading condition 
applied. The geoaetry of the lap-shear joint studied In this section Is 
Identical to that In the previous cases. The two adherends are the same 
aluminum, and the adhesive Is an epoxy resin with material constants 
reported previously. Solutions obtained for various crack lengths are 
reported In Figs. 11 and 12, In which K and G are related to a/L. Both 
Kj and shown in Fig. 11 Increase rapidly for shorter cracks (e.g., 

a/L < 0.02), and change almost linearly as a/L becomes longer. The 
opening-mode stress intensity factor is approximately four times higher 
than that of the shearing-mode In the entire range of a/L studied. Effects 
of the crack lengt^h on energy release races are shown in fig. 12. Well 
defined relationships are observable. In the case of a very short crack, 
i.e., a 0 , there exists a strong Interaction between the crack tip and 
the edges of the joint, Basic nature of the short interfacial crack and 
its interaction with geometric boundaries are currently under Investigation. 
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5. SUMMARY AND CONCLUSIONS 


An investigation of the interfacial crack behavior in adhesively 
bonded lap-shear joints is presented. The method of analysis is formu- 
lated on the basis of conservation laws in elasticity of nonhomogeneous 
solids and fundamental relationships in fracture mechanics of interface 
cracks. The current approach provides a convenient and accurate means to 
examine the basic nature of Interface cracks in adhesive joints. Fracture 
parameters such as stress intensity factors and associated energy release 
rates describing crack- tip deformation and stresses are determined. Complex 
failure modes in the joints are studied. Solutions are obtained for problems 
with various kinds of adhesive and adherends, joint configurations and 
interfacial crack lengths. Based on the results obtained, the following 
c one lus ions may be drawn; 

1. Adhesive properties have significant effects on the Interfacial 

crack tip response. In the lap-shear joint with a given adherend, the 
opening-mode stress intensity factor is always larger than that of the 
shearing-mode. The difference between K^. and increases with 

E^/E 2 . ^^hile both and Kj,^ decrease with decreasing E 2 , the total 

energy release rate and that associated with the opening-mode fracture 
increase rapidly - a phenomenon unique to the interface crack and not 
observed in homogeneous crack problems in general. 

2. The change of two identical adherends to dissimilar ones influences 
appreciably the failure behavior of the joints. Kj. increases rapidly but 

K _ remains relatively constant as the stiffness of the upper adherend 
approaches that of the lower one. An opposite trend in the changes of 
Kj and is observed as > 20. As Ej^/E 2 exceeds 20, crack 
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ilstance of the Joints hecooies relatively insensitive to the dissimilarity 
the two adherends. The Joint has a higher fracture resistance as the 
ler adherend becomes lesiit stiff than the lower one. 

3. Effects of the adhesive layer thickness on interfacial fracture are 

iferent from those on a center-of-the-bond cohesive crack in lap-shear 
:nts. is shown to be insensitive to the change of the adhesive thick- 

is, but Increases appreciably. The decrease of the adhesive thick- 

ts leads to an Increase of fracture resistance in the Joint due to reduc- 
m of crack extension driving force. 

4. Increasing the interfacial crack length Increases and 
Dultaneously . However, stress intensity solutions increase more rapidly 
: shorter cracks than for longer ones. For the cases of a/L > 0.02, 
)roximately linear relationships between K's and a/L are observed. 

a semi- logarithmic plot, Gj , G^j. and G vary with a/L almost linearly 
), suggesting that well defined relations among them may exist. 
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Fig. 4 Gaussian Stations, Finite Element Mesh and Path for M Integr^ Around Interface Crack Txp 
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Fig. 5 Stress Intensity Factors for Interface Crack in Lap Shear Joints with 
Various Adhesive Moduli (E.= E = lOxlO® psi, v = v = 0.33) 



Energy Release Rates for Interface Crack in Lap Shear Joints with 
Various Adhesive Moduli (E-= E,* 10>cl0® psi, v-» 0.33) 
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Fig. 8 Effects of Adherend Properties on Energy Release Rates for Interface Crack In Lap Shear Joints 
(E = 0.5x10^ psi, v,»0.35, E^* 10x10® psi, 0.33) 



Effects of Adhesive Layer Thickness t 2 on Crack-Tip Stress Intensity Factors in Lap Shear 
Joints (ti= to= 0.05 in, E,= 10x10® psi, v^= Vo= 0.33, E 2 = 0.5x10® psi, V 2 = 0.35) 
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Fig. 10 Effects of Adhesive Layer Thickness t2 on Energy Release Rates 
In Lap Shear Joints (tj^» t3* 0.05 In, Ei* £3* 10x10® psl, 

Vi* ^ 3 “ 0.33, E 2 " 0 . 5 x 10 ® psl, V 2 * 0.35) 




Fig. 11 Stress Intensity Factor Solutions Vs. Interface Crack. Length 
in Lap Shear Adhesive Joints (L » 0.5 in., E 3 « 10x10® psi, 
Vi“ V 3 “ 0.33, E 2 * 0.5x10® psi, V 2 “ 0.35) 
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Fig. 12 Energy Release Rates Vs. Interface Crack Length in Lap Shear 
Adhesive Joints (L * 0,5 in., Ei« E 3 * 10x10® psi, 

Vi« V 3 - 0.33, E 2 “ 0.5x 10® psi, V 2 * 0.35) 




